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Abstract. Ageing phenomena are observed in a large variety of dynamical systems 
exhibiting a slow relaxation from a non-equilibrium initial state. Ageing can be 
characterised in terms of the linear response R{t, s) at time t to a local perturbation 
at time s < t. Usually one distinguishes two dynamical regimes, namely, the quasi- 
stationary regime, where the response is translationally invariant in time, and the 
ageing regime, where this invariance is broken. In general these two regimes are 
separate in the sense that the two limits of (a) taking s — > cxd while keeping t/s 
fixed, and (b) taking i, s — > oo with fixed t — s, give different results. In recent years, 
ageing was also investigated in the context of homogeneous critical systems such as 
the Glauber-Ising model and the contact process. Here we argue that, in contrast 
to a widespread believe, homogeneous critical systems do not have a separate quasi- 
stationary regime, hence the two limits do commute. Moreover, it is discussed under 
which conditions two particular exponents, denoted as a and a' in the literature, have 
to be identical. 
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1. Introduction 

In statistical physics a dynamical system prepared in a certain initial state is said 
to exhibit ageing if its physical properties vary slowly with time [1]. Such a time 
dependence is typically observed in disordered systems such as spin glasses which 
memorise their initial state for very long time. As a hallmark of ageing, one observes 
that quantities such as two-point correlation functions are not translationally invariant 
in time. 



Ageing can be characterised in various ways, e.g. by studying the autoresponse 
function R{t, s). This function quantifies how a weak perturbation by an external field h 
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Figure 1. Cartoon of the autoresponse function R{t, s) measured in a 
disordered system (see text). 



at time s exerted locally at position x will change the averaged order parameter s)) 
at some later time t > s: 



Fig. [T] sketches how the autoresponse function is measured. At t = the system is 
quenched to a non-equilibrium initial state. As time evolves the system relaxes slowly 
towards the stationary state (without being able to reach it) , accompanied by a gradual 
change of the averaged order parameter, as indicated by the changing background colour 
in the figure. After a certain waiting time s a local external field is switched on for a 
short time, imposing a local change of the order parameter (e.g. a spin flip). As 
time proceeds this perturbation spreads in space and becomes less intense, as expressed 
schematically by the parabola-shaped region in the figure. The autoresponse function 
is just the average variation of the order parameter measured at time t > s at the same 
position from where the perturbation originated. 

In numerous studies of disordered systems it was observed that the temporal 
evolution after the perturbation can be divided into two different dynamical regimes 
(see e.g. Zippold et al. [2]). Initially the influence of the external field dominates and 
causes a response which is essentially insensitive to the underlying slow relaxation of 
the entire system. In this so-called quasi- stationary regime the response function is 
translationally invariant in time, i.e. it depends to leading order ont — s but not on s: 



After some time, however, the response becomes weaker and begins to 'feel' the 
relaxation of the system. In this so-called ageing regime the response function Rag{t, s) 
is no longer translational invariant in time, rather it depends on both parameters t and 
s in a non-trivial way. Both regimes are connected by a smooth crossover taking place 




(1) 




(2) 
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at some typical time scale Tc(s) > s. Surprisingly, in disordered systems this crossover 
time was found to scale as [2] 

Tc{s) - s ~ s'^ (3) 

with an exponent < 7 < 1. In other words, the temporal size of the quasi- stationary 
regime grows with the waiting time s but it grows slower than s itself. 

As noted by Zippold et ai, Eq. ([3]) has the non-trivial consequence that the two limits 
of 

(a) first taking t, s — > 00 with C, = t/s fixed and then sending ^ ^ 1, and 

(b) sending t,s —>■ 00 while keeping the difference t — s fixed 

are generally different. This inequivalence is illustrated in Fig. [21 where the two 
dynamical regimes are sketched in a double-logarithmic t-s-plane. The typical path 
according to limit (a) is marked as a green line. As can be seen, this line runs parallel 
to the diagonal and enters the ageing regime after some time. When s, t are taken 
to infinity (moving along the diagonal) the width of the quasi-stationary regime in 
the log- log representation shrinks to zero. Therefore, the second step of (a), the limit 
t/s — > 1, takes place entirely within the ageing regime and hence this limit is determined 
exclusively by the properties of the ageing regime. Contrarily, in the case of limit (b), 
where the difference t — s is kept constant, the path (red line) eventually enters the 
quasi- stationary region, hence this limit will be determined by the properties of the 
quasi- stationary regime only. For this reason it is plausible that the two limits give 
generally different results. 

The interest in ageing phenomena stems from the fact that the ageing regime 
displays scale-invariant features. This means that a change of time scales s —>■ bs, t ^ bt 
by a factor 6 > can be compensated by an appropriate rescaling of the response 
function R b^^^"'R, where a is a certain exponent depending on the scaling dimension 
of the order parameter and the response field. Scale invariance in the ageing regime 
implies that the response function obeys the scaling form 

i?ag(t, S) ~ S-'-'^fRit/s) , t > T,{S) (4) 

where fn is a scaling function. For large arguments it is found that the scaling function 
decays algebraically as 

(5) 

where Xr is a new non-trivial exponent, the so-called ageing exponent, while z denotes 
the dynamical exponent of the ageing regime. The scaling function and the exponents 
a,XR,z are expected to be universal, i.e., they are determined by the symmetry 
properties of the model irrespective of its microscopic details. 
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Figure 2. Inequivalence of tiie limits (a) and (b). According to [2] the response 
to a perturbation can be split into a quasi-stationary regime and an ageing regime, 
connected by a crossover at a typical time scale Tc{s) > s (dashed blue line), where 
Tc{s) — s ~ s''' with < 7 < 1 (here 7 = 1/2). The green and the red line symbolise 
typical paths along which the limits (a) and (b) are taken (see text). 



In recent years the concept of ageing was also applied to homogeneous critical 
systems with a continuous phase transition [3] such as the 2D Ising model [4] and 
directed percolation [5-7] without quenched disorder. However, these models differ 
from previously studied systems such as spin glasses in so far as they have a much 
simpler scaling behaviour. The purpose of this contribution is to point out that in this 
case - contrary to a widespread believe - the two dynamical regimes are not separate 
but connected by an ordinary scale-invariant crossover, meaning that the two limits (a) 
and (b) give identical results. 

The second issue addressed in this paper concerns a recently proposed theory of 
local scale invariance (LSI), which introduces an additional exponent a'. The arguments 
presented here imply that for the 2d Ising model and the contact process one has a = a'. 

2. Ageing in systems at criticality 

In what follows let us consider a stochastic process on an infinite lattice with the 
following properties: 

(i) The model is homogeneous in space and time (no quenched disorder) and exhibits 
a continuous phase transition. 

(ii) The transition is characterised by simple scaling (as opposed to multiscaling) 
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involving three critical exponents, namely, the order parameter exponent /3, the 
correlation length exponent z/ = z/^, and the dynamical exponent z = vw/vx- 

(iii) The system has an autonomous response in the stationary state, i.e., the response 
function R{t, s) at criticality is non-zero in the limit s ^ oo with t — s fixed. 

Examples of such models are the critical Ising model in two dimensions (2d) with 
standard heat bath dynamics as well as critical directed percolation in 1+1 dimensions. 
These models display asymptotic scaling laws which are well understoodl^l Clearly, these 
scaling laws are valid in the stationary state at criticality and thus they will determine 
the scaling properties in the quasi-stationary regime. However, in contrast to disordered 
systems such as spin glasses, the same scaling laws are also valid in the ageing regime 
t ^ s. As will be shown below, this implies that the scenario suggested by Zippold et 
al. [2] does not hold for such systems. 

To see this let us assume the existence of two dynamical regimes separated by a 
crossover at a typical time scale Tc{s) — s s'^ with 7 < 1. In the quasi-stationary 
regime, where the response function Rst{t,s) = Rst{t — s) is translational invariant in 
time, invariance under rescaling 

t^Xt, s-*Xs, x^X^/'x, R^X-f^/^'R (6) 

implies that R{Xt — As) = X~^^'^^R{t — s), leading to the well-known result 

i?,t(t, s) ~ ro (t - s)-'^/^^ . t « r,(s). (7) 

Note that for models with an autonomous response in the stationary state, the 
proportionality constant tq does not depend on s. 

In the ageing regime, the system is still invariant under the same type of scale 
transformation ([6]) but translational invariance in time is broken, leading to the weaker 
scaling form 

i?,g(t,s) ~ s-''^'''fR{t/s), t»r,(s). (8) 

hence a = P/uz — 1 in Eq. (jll). This scaling form has been verified for various systems 
and is undisputed in the literature. 

Both scaling forms are smoothly connected at the crossover t = Tc(s), that is 
Rst{t,Tc{s)) ^ i?ag(t, Tc(s)), leading to the condition 

rofniQ ^ {^,-ir^/''\ (9) 

where '■= Tc{s)/s. For 7 < 1 the argument C,c varies with s, fixing the functional 
form of the scaling function fn. However, this particular functional form would imply 

t Note that any lattice model displays corrections to scaling caused by lattice effects on short distances 
and time scales (denoted as Tmicro in the literature). Throughout this paper we assume that distances 
and time scales are large enough so that such lattice effects can be neglected. 
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translational invariance Rag{t,s) = Rag{t — s) in parts of the ageing regime, which is 
definitely not true. Therefore, the scenario suggested by Zippold et al. does not hold for 
critical systems such as the 2d Ising model or directed percolation. The only way out is 
to set 7 = 1 so that Eq. ([9]) determines the scaling function only in a single point. This 
means that the crossover time Tc(s) ~ s scales like any other time scale, which is quite 
natural in the context of second-order phase transitions. In the representation of Fig. [2] 
the ageing regime would be a strip of constant width. As can be seen, this immediately 
implies that the two limits (a) and (b) will give identical results. 

3. Test of Local Scale Invariance 

Recently a theory of Local Scale Invariance (LSI) [3, 8] has been applied to various 
homogeneous many-particle systems with second-order phase transitions. This theory 
assumes that the linear response function formed from so-called quasi-primary fields 
transforms covariantly under the action of a group of local scale transformations. Like 
conformal invariance, this extended symmetry determines the functional form of scaling 
functions. In particular, LSI predicts that the scaling function /r(^) factorises into two 
power laws (see Eq. (26) of Ref. [8]) 

/k(0 = /oei-'-W^(e-l)-i-'. (10) 

with an additional exponents a' which may be different from a. However, knowing that 
the aforementioned limits commute, this scaling form would imply that for t/s —>■ 1^ 
the response function behaves as 

R{t,s) ~ /o(ts'^'-'^-s)"'"'' . (11) 

Obviously, this behaviour is in contradiction with Eq. ([2]) unless a' = a or fo = 0. 

Originally the theory of LSI was formulated for systems with translational 
invariance in time, which led to the constraint a' = a. The initial hope was that 
this theory is fundamental enough to describe a broad range of apparently unrelated 
models. However, later it became clear that the critical Ising model in 2d and the 1+1- 
dimensional contact process do not obey the predicted scaling forms for the autoresponse 
function [9-11]. As a way out, Henkel and collaborators generalised the theory, relaxing 
the symmetry by giving up translational invariance in time which allows a and a' to 
be different. Fitting numerical data of various critical models to the predicted response 
function they found a good agreement when [8] 

(—0.187(20) for the 2d Ising model with heat bath dynamics, 
—0.022(5) for the 3d Ising model with heat bath dynamics, 
+0.270(10) for the 1+1-dimensional contact process. 

However, these models have an autonomous response (/o > 0), hence the above 
arguments require that a' = a. This suggests that the numerical estimates reported 
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above have no deeper physical meaning; they merely result from the attempt to fit a 
wrong scaling function to numerical data within a limited window of ^. 

The proponents of generahsed LSI refer to the circumstance that for the Id Ising 
model one can prove that a' — a = —1/2. However, this model becomes critical at zero 
temperature with a non-fluctuating steady state so that it does not have an autonomous 
response in the limit s ^ oo, hence /o = 0. Therefore, it is no surprise that a ^ a' . 

4. Conclusions 

Ageing was originally introduced as a concept to describe the propagation of pertur- 
bations in disordered systems. As a key property, such systems show two separate 
dynamical regimes, namely, a quasi-stationary one and the so-called ageing regime. 

In homogeneous critical systems such as the 2d Ising model and the contact 
process, the situation is very different. Here the entire dynamical range (apart from 
corrections due to lattice effects) is determined by the same type of scaling. In fact, 
in homogeneous critical systems a quenched initial state plays the role of a temporal 
boundary condition analogous to a surface in a magnetic system at criticality. Such a 
boundary condition induces one additional surface exponent (here Xr) but it neither 
generates additional scales nor pushes the system into a different universality class. 
Instead of two dynamically different regimes homogeneous critical systems display a 
simple crossover phenomenon which can be described in the framework of standard 
scaling theory. 
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